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We have measured the neutron structure function g2 and the 
virtual photon-nucleon asymmetry A?> over the kinematic range 
0.014 < x < 0.7 and 1.0 < Q 2 < 17.0 by scattering 48.3 GeV lon- 
gitudinally polarized electrons from polarized 3 He. Results for A?> 
are significantly smaller than the \fR positivity limit over most of 
the measured range and data for gJ> are generally consistent with 
the twist-2 Wandzura-Wilczek prediction. Using our measured g?; 
we obtain results for the twist-3 reduced matrix element d%, and 
the integral $%£{x)dx in the range 0.014 < x < 1.0. Data from this 
experiment are combined with existing data for g^ to obtain an 
average for dV^ and the integral f^{x)dx. 



1 Permanent Address: Institut des Sciences Nucleaires, IN2P3/CNRS, 38026 Greno- 
ble Cedex, France 

2 Permanent Address: University of Bochum, D-44780 Bochum, Germany 

3 Present Address: Florida International University, Miami, FL 33199 



2 



The deep inelastic spin structure functions gi(x, Q 2 ) and g2^x,Q 2 ), which de- 
pend on the Bjorken scaling variable x and the virtual photon four-momentum 
squared — Q 2 , provide insight into the internal spin structure of the nucleon. 
A large set of data for g x now exists for the proton, deuteron [1,2] and neu- 
tron [3,4]. These data have been used to test the fundamental Bjorken sum 
rule, and within the framework of the quark-parton model (QPM), to mea- 
sure the quark contribution to the nucleon's spin. The g2 structure function 
contains contributions from both the longitudinal and transverse polarization 
distributions within the nucleon. It is sensitive to higher twist effects such 
as quark-gluon correlations and quark mass contributions, and is not easily 
interpreted in the QPM where such effects are not included. However, by in- 
terpreting g2 using the operator product expansion (OPE) within QCD [5,6], 
it is possible to study contributions to the nucleon spin structure beyond the 
simple QPM. 

The OPE allows us to write the hadronic matrix element in deep inelastic 
scattering (DIS) in terms of a series of renormalized operators of increasing 
twist [5,6]. The leading contribution is twist-2, with higher twist terms sup- 
pressed by powers of l/Q. Keeping only terms up to twist-3, the moments of 
gi and g 2 at fixed Q 2 can be related to the twist-2 and twist-3 reduced matrix 
elements, aj and dj [6], 

i 

J x^g 1 (x,Q 2 )dx = ^-, j = 0,2,4,... 
o 

Jx i g 2 (x,Q 2 )dx = ±^—(d j -a j ), j = 2,4,... (1) 
o ^ 

In the expressions above, dj directly appears in the equation for g2 allowing 
us to study the higher twist structure of the nucleon at leading order. An 
expression for the twist-2 part of g 2 was derived by Wandzura and Wilczek [7] 
from these sum rules assuming that the twist-3 contributions dj, are negligible, 

g? w (x, Q 2 ) = - gl (x, Q 2 ) + / ^-l^W. (2) 

j x 

x 

Comparing measured values of g 2 with this prediction enables us to extract 
information about higher twist contributions to g 2 . There is an additional 
twist-2 contribution to g 2 [8,9] beyond the g\ vw term which arises from the 
transverse polarization density in the nucleon, Ht{x, Q 2 )- However, this term 
is suppressed by the ratio of the quark to nucleon mass m/M in DIS [8] and 
will be neglected in this analysis. 
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The structure function g2(x, Q 2 ) may be expressed in terms of two measur- 
able asymmetries, -Am (a;, Q 2 ) and A±(x, Q 2 ), corresponding to longitudinal and 
transverse target polarization with respect to the incoming electron beam he- 
licity, 



g2(x,Q 2 ) 



F 2 (x,Q 2 )(l + 7 2 



V 



2x [1 + R(x, Q 2 )} 2d sin # 



A 



E + E'cos 9 
E> 



Aw sintf 



, (3) 



where E and E' are the incident and scattered electron energies, is the 
scattering angle, 7 = 2Mx/y/Q^, y = (E - E')/E, d= (1 - e)(2 - y)/y[l + 
eR(x, Q 2 )], and e" 1 = 1 + 2 [1 + 7~ 2 ] tan 2 (#/2). Fits to existing data were used 
for the unpolarized structure function F 2 (x, Q 2 ) [10] and for R(x, Q 2 ) [11], the 
ratio of longitudinal to transverse virtual photon absorption cross sections. 
At small scattering angles, the term An sin# is small, and consequently the 
dominant contribution to g2 comes from A±. 



Spin dependent DIS can also be described in terms of the spin asymme- 
tries Ai(x,Q 2 ) and A 2 (x, Q 2 ) for virtual photon absorption. The asymmetry 

A 2 (x, Q 2 ) is bounded by the positivity limit |A 2 (x, Q 2 )\ < ^jR(x, Q 2 ), and like 
g2, it is dominated by A±, 



A 2 (x,Q 2 ) = 



7(2 -y) 
2d sin 6 



■ y(l + xM/E) 

A± (1-1/) +A " 



sin 9 



(4) 



Measurements of g 2 and A 2 exist for the proton [12,13] and deuteron [12], and 
in the case of the neutron, a measurement was made at the Stanford Linear 
Accelerator Center (SLAC) by scattering 26 GeV polarized electrons from 
polarized 3 He [3]. In this Letter, we report a new measurement of g 2 and A 2 
for the neutron made during experiment E154 at SLAC. For this experiment, 
the beam energy was increased to 48.3 GeV and two new large-acceptance 
spectrometers were constructed to provide broader kinematic coverage than 
previously measured. Results from this experiment for g™ and A™ have been 
reported elsewhere [4], and we focus here on the measurement of A± and the 
subsequent determination of g 2 and A^. 

The target was a 30 cm long, thin-walled glass cell containing approximately 
10.5 atmospheres (as measured at 20°C) of 3 He gas. The helium nuclei were 
polarized by spin-exchange collisions with rubidium atoms that were polar- 
ized by optical pumping [14]. The system was designed to allow continuous 
pumping of the target polarization in the longitudinal direction only. There- 
fore, to obtain transverse polarization, the 3 He spins were first pumped to a 
longitudinal polarization of 48% and then rotated to the transverse direction 
using two orthogonal sets of Helmholtz coils. In the transverse orientation, the 
polarization decreased to 33% over a period of 24 hours, at which time the 
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target was re-polarized. Approximately 7 x 10 6 electron events were recorded 
during two cycles of transverse running. 



The electron beam was produced in 250 ns pulses at a rate of 120 Hz, each con- 
taining approximately 3 x 10 10 electrons. The average beam polarization was 
measured to be 0.826 ±0.023 using a M0ller polarimeter [15], and the helicity 
of each pulse was chosen randomly to reduce helicity-dependent systematic 
errors. Scattered electrons were measured using two independent spectrome- 
ters at scattering angles centered around 2.75° and 5.5°, and the asymmetry 
Aj_ was calculated as 



where N + and N~ are the measured electron rates for positive and negative 
beam helicities corrected for detection efficiency and normalized to the incident 
charge, P^ and P t are the beam and target polarizations, and / is the dilution 
factor which corrects for electrons that scattered from materials in the target 
system other than 3 He. 

Our measured asymmetry included not only DIS events, but also pions mis- 
identified as electrons, and electrons produced in charge symmetric hadron 
decays. The rates and asymmetries for these backgrounds were measured and 
used to correct A±. The asymmetry was also corrected for internal [16] and 
external [17] radiative effects. Uncertainties in the radiative corrections were 
estimated by varying the input models over a range consistent with the mea- 
sured data. Finally, a neutron result was extracted from Aj_ by applying a 
correction for the 3 He nuclear wave function [18] and using g^ w obtained 
from a fit to existing proton data [1,2] for g{. 

Results for A2 and gJJ from both spectrometers are given in Table 1 with sta- 
tistical and systematic errors. The data cover the kinematic range 0.014 < 
x < 0.7 and 1.0 < Q 2 < 17.0 (GeV/c) 2 with an average Q 2 of 3.6 (GeV/c) 2 . 
Systematic errors are dominated by the uncertainty in the radiative correc- 
tions, but are significantly smaller than the statistical error over the entire 
data range. No evidence of Q 2 dependence was seen for A2 or g£ within the 
experimental errors and the data from both spectrometers were averaged. The 
results for Ag are shown in Fig. 1 along with the y/R positivity limit and pre- 
vious data from SLAC experiment E142 [3]. The data are in good agreement 
with the E142 measurement and are significantly smaller than the positivity 
limit over most of the measured range. Results for xg% are shown in Fig. 2 
along with the twist- 2 prediction, xg^ w . To calculate g^ w , we assume that 
gi/F\ is independent of Q 2 , and use a fit to our measured g™ data [4]. A com- 
parision of our data with g^ w over the measured range gives a xV(dof) of 
1.02 indicating good overall agreement. However, the data clearly do not rule 
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Table 1 

Results for A?> and g 2 for the 2.75° and 5.5° spectrometers. 



x range 


< x > 


<Q 2 > 


An 










(GeV/c) 2 


±stat±syst 


±stat±syst 








2.75° 


Spectrometer 




0.014 - 


0.02 


0.017 


1.2 


0.03 ±0.07 ±0.01 


7.36 ±15.74 ±2.24 


0.02 - 


0.03 


0.025 


1.6 


0.00 ±0.06 ±0.01 


0.15 ±7.19 ±0.98 


0.03 - 


0.04 


0.035 


2.1 


-0.11 ±0.06 ±0.01 


-7.90 ±4.91 ±0.96 


0.04 - 


0.06 


0.049 


2.6 


0.10 ±0.06 ±0.01 


4.60 ± 2.50 ± 0.54 


0.06 - 


0.10 


0.077 


3.4 


0.06 ±0.06 ±0.01 


1.32 ± 1.34 ±0.25 


0.10 - 


0.15 


0.122 


4.1 


0.13 ±0.11 ±0.03 


1.22 ±0.95 ±0.24 


0.15 - 


0.20 


0.173 


4.7 


-0.03 ±0.18 ±0.03 


-0.08 ±0.81 ±0.14 


0.20 - 


0.30 


0.242 


5.1 


-0.25 ±0.24 ±0.05 


-0.48 ±0.51 ±0.11 


0.30 - 


0.40 


0.341 


5.6 


0.63 ±0.55 ±0.13 


0.54 ±0.46 ±0.15 


0.40 - 


0.50 


0.425 


5.9 


0.16 ± 1.40 ±0.04 


0.04 ±0.57 ±0.02 








5.5° 


Spectrometer 




0.06 - 


0.10 


0.084 


5.5 


0.16 ±0.10 ±0.02 


4.08 ± 2.40 ± 0.43 


0.10 - 


0.15 


0.123 


7.2 


0.01 ±0.08 ±0.02 


0.23 ± 1.00 ±0.20 


0.15 - 


0.20 


0.173 


8.9 


0.05 ± 0.11 ± 0.02 


0.40 ± 0.72 ± 0.15 


0.20 - 


0.30 


0.242 


10.7 


0.15 ±0.14 ±0.03 


0.48 ±0.41 ±0.10 


0.30 - 


0.40 


0.342 


12.5 


-0.21 ±0.27 ±0.03 


-0.22 ±0.31 ±0.04 


0.40 - 


0.50 


0.442 


13.8 


-0.36 ±0.53 ±0.05 


-0.16 ±0.24 ±0.03 


0.50 - 


0.70 


0.564 


15.0 


-0.04 ± 0.96 ± 0.06 


-0.01 ±0.13 ±0.01 



out the possibility of large twist-3 contributions and show marginal agreement 
with the twist- 2 prediction in the region 0.03 < x < 0.1. 

In order to quantify the possible contribution of higher twist effects to gg, 
Eq. 1 can be solved for the twist-3 reduced matrix elements (X] at fixed Q 2 , 



i 

d](Q 2 ) =2jx 3 



dx, j = 2,4,... (6) 



The combination of gi and g2 in the above expression effectively cancels any 
twist-2 components allowing us to look for a net twist-3 contribution to g 2 . 
The matrix element was calculated using our g 2 (x,Q 2 ) data and a fit to our 
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Fig. 1. The asymmetry A?> for both spectrometers combined, and the corresponding 
\fR positivity limits. Also shown are the neutron data from SLAC experiment E142. 
Errors are statistical only 
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Fig. 2. The structure function xgJf for both spectrometers combined. Also shown is 
the the twist-2 g^ W prediction. Errors are statistical only. 

measured gi(x, Q 2 ). Because the integrand in Eq. 6 is purely twist-3, we as- 
sumed the unmeasured region 0.7 < x < 1 behaves like (1 — x) 2 as suggested 
by Brodsky et al. [19], and fit our last data point to this form to extrapolate to 
x — 1. We neglected any contribution from the region < x < 0.014 because 
it is suppressed by the x^ term. For the j = 2 moment, we obtain a value of 
dl = -0.004 ± 0.038 ± 0.005 with an average Q 2 of 3.6 (GeV/c) 2 . The contri- 
bution from the high-x extrapolation is much smaller than the experimental 
errors and does not significantly change the result for the matrix element. 

Data from SLAC experiments E142 [3] and E143 [12] were combined with this 
experiment to yield a average neutron result for g 2 with an average Q 2 of 3.0 
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• Combined SLAC data 

- Xg 2 ^w 

' 1,1,, 
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> 



Fig. 3. The structure function xg% for SLAC experiments E142, E143 and E154 com- 
bined. Also shown is the twist-2 gY W prediction. The average Q 2 for the combined 
data is 3.0 (GeV/c) 2 . Errors are statistical only. 

Table 2 

Comparison of experimental and theoretical results for the reduced twist-3 matrix 
element d%- 





dl x 10 2 


Q 2 (GeV/c) 2 


El 54 result 


-0.4 ±3.8 


3.6 


SLAC Average 


-1.0 ± 1.5 


3.0 


Bag model [8] 


-0.253 


5.0 


Bag model [12,20] 


0.03 


5.0 


Bag model [21] 





5.0 


QCD sum rule [22] 


-3± 1 


1.0 


QCD sum rule [23] 


-2.7 ± 1.2 


1.0 


Lattice QCD [24] 


-0.39 ±0.27 


4.0 



(GeV/c) 2 . Neutron results were extracted from the E143 proton and deuteron 
data assuming a 5% D-state in the deuteron. The results are shown in Fig. 3 
along with the g^ w prediction. Comparing the combined data with g^ w 
gives a x 2 /(dof) of 1.01, again indicating good agreement with g^ w ■ Using 
the combined data, we obtain the result d% = —0.010 ± 0.015 at an average 
Q 2 of 3.0 (GeV/c) 2 . The measured d% along with various model predictions 
are summarized in Table 2, and while the data are consistent with zero, the 
precision is insufficient to rule out models which contain significant twist-3 
contributions. 

From the OPE it is not possible to obtain an expression for the j — moment 
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of g2- However, Burkhardt and Cottingham [25] have derived the sum rule 
Jq 1 g2(x)dx = 0, which is valid to first order in pQCD [26], using dispersion 
relations for virtual Compton scattering. To evaluate the integral, the g^ w 
expression in Eq. 2 was used to evolve the twist-2 part of our measured g 2 
to a Q 2 of 3.6 (GeV/c) 2 assuming gi/iq is independent of Q 2 and fitting our 
g" data as before. At large x, we see from Eq. 2 that gY w ~ — gi and we 
therefore assume that g2 oc (1 — x) 3 for the extrapolation to x — 1. The result 
is Jo m4 g 2 (x)dx = 0.19 ± 0.17 ± 0.02 with an average Q 2 of 3.6 (GeV/c) 2 . 
The Q 2 evolution and high-x extrapolation do not contribute significantly 
to the integral and the uncertainties in these quantites are included in the 
error. Combining this result with data from SLAC experiments E142 [3] and 
E143 [12] yields a result of / . 014 g 2 (a;)dx = 0.06 ± 0.15 at an average Q 2 of 
3.0 (GeV/c) 2 , which is consistent with the Burkhardt-Cottingham sum rule. 
However, this does not represent a conclusive test of the sum rule because the 
behavior of g% as x — > is not known. 

In summary, we have presented a new measurement of and g£ in the 
kinematic range 0.014 < x < 0.7 and 1.0 < Q 2 < 17.0 (GeV/c) 2 . Our results 
for Ag are significantly smaller than the \fR positivity limit over most of the 
measured range and data for gJ? are generally consistent with the twist-2 g^ w 
prediction. The values obtained for the twist-3 matrix element d% from this 
measurement and the SLAC average are also consistent with zero. However, 
further measurements are needed to make a conclusive statement about the 
higher twist content of the nucleon. 
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